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Perovskite solar cell (PSC) is one of the advanced third-generation solar 
cells that have rapid efficiencies but instability issues in terms of air, 
moisture, and UV light sensitivity becomes a barrier to commercialization. 
The instability issue is due to the charge accumulation at the interface of the 
PSC reducing its efficiency. This research focuses on the operation of PSCs 
through a one-dimensional (1D) drift-diffusion model for planar 
heterojunction PSCs. The model also accounts for the electric potential by 
Poisson’s equation, the ion generation, and the recombination rate. The 
method of lines technique is applied to solve the model for the perovskite 
layer numerically using the finite difference method which is then solved 
forward in time using the ‘ode/5s’ solver in MATLAB. It is highlighted that 
the comparison with the experimental data from the reference shows good 
agreement. The effect of parameter thickness variation of the perovskite 
layer upon the efficiency of PSCs is analyzed. The result shows that the best 
efficiency obtained is 19.77% obtained at thickness 0.25 um. The results 
may prove useful for a guideline of the cell thickness that predicts the 
perovskites cell performance. 
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1. INTRODUCTION 


The solar cell is a device that converts light energy into photovoltaic energy. The usage of solar cells 
that are pollution-free, low-cost, and renewable energies are the main advantages that contribute to a positive 
effect on living and economics [1], [2]. Perovskite solar cell (PSC) is one of the advanced technologies 
developed from the dye-sensitized solar cell (DSSC) that is categorized as a third-generation solar cell [3], 
[4]. Since 2009, the perovskites solar cell has been remarkable due to cost-effectiveness, high-power 
conversion efficiency, strong optical absorption, and ease of fabrication [5]. 

However, the stability issue in terms of moisture, light, toxicity level, and thermal lead to current- 
voltage hysteresis, reducing the efficiency of the PSC [6]-[8]. The instability issue is related to the 
mechanisms of charge transport via diffusion (charge density gradient) and advection (by an electric field) 
and charge release at the interface between layers during ion migration which dominates the global electrical 
field [8]-[10]. The movement of the charge transport is related to the dynamic physics of the PSC which 
requires a proper mathematical model that describes charge transport layers [11]. 
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Various works perform the modeling of the charge accumulation at the interface of the PSC. Based 
on the finding from Neukom et al. [12] considering the ion and charge traps that are capable of relating the 
hysteresis behavior in PSCs using experimental and numerical simulation Fluxim but the mobility of ion 
migration is assumed to be constant and the concentration of upper limit boundaries does not being imposed. 
The impedance spectroscopy technique is applied by Khan et al. [13] to probe the interfacial surface of the 
perovskite layer to study the interface charge accumulation, recombination losses, and charge transport 
mechanism but the technique is highly sensitive. The finding shows that the organohalide type MAPBI3 
contributes to low interfacial recombination compared to organohalide ((FAPbI;)o.g5 (MAPDBr3)945) [13]. 
The equivalent electrical circuit for modeling the PSC involves the series of two diodes that represent the 
active junction at the holes transport layer and electrons transport layer as the perovskite layer absorbs the 
sunlight energy [14], [15]. Using equivalent circuit modeling to study the physics of PSC is too 
straightforward and lacks the electrochemical behavior of PSC while the atomistic density functional theorem 
(DFT) requires high computational cost and applies to a limited number of atoms [16]. The work from 
Courtier et al. [16] shows that the transport layer properties affect the performance of the PSC depending on 
the degree of the hysteresis using the extension of the numerical scheme from. To fill in the gap, the drift- 
diffusion model is found to be the best model for modeling the PSC due to the ability to account for the 
dynamic physics of the charge transport model involving the drift and diffusion mechanisms and potential 
difference with extremely small Debye layer [16], [17]. 

The solution to the drift-diffusion equation may vary such as using a solar cell simulator 
(SCAPS-1D [18], [19], wxAMPS [20], IonMonger [21]) directly or using the numerical method and 
asymptotic analysis [22], [23]. The challenge that arises when using simulators is the generation profile 
needed using the transfer matrix method [19]. The drift-diffusion equations consisted of the coupled ordinary 
differential equations (ODE) and partial differential equation (PDE) [22], [24]. Instead of using the complex 
asymptotic analysis [16], [17], the method of line (MOL) technique is proposed to solve the drift-diffusion 
equations which is also the main contribution in this research. The method of the line has advantages in 
decoupling the spatial and temporal discretization processes [25]. The result is compared with the reference 
[24] for validation purposes. Next, the thickness of the perovskite layer (absorber) is varied to analyze the 
effect on the efficiency of the PSC. The main contribution to the research is the MOL technique to solve the 
perovskites solar cells model. 


2. MATHEMATICAL MODEL AND NUMERICAL PROCEDURE 
This section discusses the mathematical model of the PSC and the numerical procedure to solve the 
model. 


2.1. Modeling of the perovskite solar cell 

Figure 1 shows the schematic diagram of solar cells in which the perovskite is sandwiched between 
the holes transport layer (HTL) and electrons transport layer (ETL). The holes transport layer (HTL) is made 
up of spiro-OMETAD and the ETL is made up of titanium oxide [24]. The p and n represent holes and 
electrons, respectively, where L represents the perovskite thickness. 


Donor phase (p-type) Blend phase Acceptor phase (n-type) 


L L L 
x= -Lli ->5 x=-3 x=5 x=lz+3 


Figure 1. The schematic diagram of the PSC. The solid line represents the highest occupied molecular orbital 


(HOMO) while the least unoccupied molecular orbital (LUMO) (dashed line). The p-type donor layer 


L L : L L L L 
(-L, = ae 5), the perovskite (- 3<%x*< 5) and the n-type acceptor layer (- <x<L,+ 5) 


The sunlight is absorbed in the perovskites layer which has a highly ordered crystalline structure 
comprised of distinct conduction band edges, Ên, and valence band edges, Pp that are divided by a band gap 
[24]. The organics materials of the donor phase (p-type) and acceptor phase (n-type) have no distinct band 
structure and are amorphous [24]. LUMO portrays the conduction band, and the energy of the acceptor is 
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denoted as $, while HOMO portrays the valence band, and the energy of the donor is denoted as $,. The 
electrons are excited into the LUMO, leaving a hole in the HOMO, resulting in the conduction process due to 
the excited electrons and hole movement between the HOMO and LUMO respectively on the adjacent 
molecules. The dynamic modeling for PSCs in this research describes the electric potential by the poisson 
equation, the charge densities and current densities for holes and electrons by drift-diffusion equations, 
continuity equations of ions, and the electron-hole pair recombination and generations are also accounted 
[17]. Therefore, the mathematical model of charged ion species (p and n) which diffuse by the ionic 


concentration gradient and advect by an electric field in the perovskite layer (-= <x< 5) in (1)-(3) [24]. 


Poisson’s equation: 
a2 
<< = 4(n—-p) (1) 


Drift-diffusion equations: 


dp , 1asP _ —n?) P=— ea 

at = q ax a(np nj X J= qDp (2 + kT zE) j 
an 1a XI, _ an _ qndg 

x go AP IT" =qDn & kT a “ 


The p,n are the density of holes and electrons, respectively, J” and J” are the current density of 
electrons and holes, ¢ is the electric potential, a(np — n?) is the recombination rate. The recombination rate 
consists of a, the radiative recombination coefficient, and, n; the doping density. The type of recombination 
considered in this work is band-to-band radiative recombination which is the recombination that naturally 
occurs in PSCs [21]. The recombination occurs when the electrons in the conduction band recombine with 
holes in the valence band resulting in the emission of a photon in between the band gap [22]. Meanwhile, q is 
the elementary charge, € is the permittivity, D, and D, are the diffusivity of electrons and holes, k is the 
Boltzmann constant and T is the temperature, respectively. The total current throughout the device is 
J=J"+ P. 

The boundary conditions at the interface of the perovskite layer to the acceptor phase, x = — =, and 


at the interface of the perovskite layer to the donor phase, x = - in (4) and (5), respectively: 


f = = E 

PiS Pv: Pla = eyg Maly 7 (4) 

A a ee _ mo 5 

Pl, at =o, n x= = Eno, Pal pat =I (5) 
where £p, En are the valence band density between the donor and perovskite layer and the conduction band 


density between the acceptor and perovskite layer, respectively. Here, 7) is an assumption of the typical 
; : PIONI L L 
charge density that carries a current, d, and @, are the potential differences at x = = and atx = A of the 


device, respectively. Here, it is assumed that the total voltage flow at the interface of perovskites at both 
phases (acceptor and donor) is the same. The potential at x = 0 is set to Ê|x=ọ = 0, which complies the 
Kirchhoff’s Voltage Law [24]. The initial condition (at t = 0) for the potential is fitted to the experimental 
data from [21]. 


B(x, 0) = Vy(661.3x7- 115.1. x° + 194.2x5 + 33.3x4-17.13x3 — 2.293x? + 0.3135x + 0.026)(6) 


The potential Vp is set equal to the thermal voltage Vy) = = which is equivalent to 20V. At x = — =, 


only the holes, p are allowed into the donor phase, which is assumed to equal to €)7o while at x = =, only the 


electrons, n exited to enter the acceptor phase and are set to Epto. 


2.2. Nondimensionalisation 
In (1) to (5) were converted into the dimensionless model and the numerical simulation was solved 
in dimensionless. The charge density, the cell thickness, and the potential are scaled to the typical charge 


carrier density, To = 7.9 x 1072m73, the width of perovskites solar cell (L), and the thermal voltage = 


respectively [21]. In this work, the hole and electron density are assumed to be the same as carrier density, TTo. 
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The superscript tilde denotes the parameter is dimensionless. The scaling form (1) to (5) is as follows (note 
that the tilde denotes the dimensionless variables): 


kT + 


ñ = non, P= Mop, E=tt, $ =F $, F= Lx, ni = 1G (7) 


Here we assume that both holes and electrons diffuse at the same diffusion rate D, = Dy = Do, 
where Dy is the typical diffusivity of the charges in which the time scaling of the cell operation is: 


rai (8) 


=p 


which leads to the following dimensionless equations: 


arg Be as 

ae aD O) 
an a an. OG ae 

or = (or nice) + yp 1) (10) 
ap _ a (ə , db pee 

Pat (2 4p) + @p-1 (11) 


a $ D3 5 
Ple: = A Pleat = Ep» fizle=zt = -1 (12) 
T br x ~ 

Plz = Vo’ alai = En Belg =1 (13) 


The dimensionless parameters in (9) to (11) are defined by: 


kT L? 
PES Y= To (14) 


Toq Do 


where À and y are the ratio of the Debye length to the width of the perovskite layer and the ratio of the typical 
charge diffusivity in bulk density to the typical diffusivity of the charges in the perovskite, respectively. The 
J-V curve is an important tool to determine the efficiency of the perovskites solar cell. The J-V curve 
function shown in (15) is used referred from reference [24], with the assumption that only the perovskite 
layer is considered. The relationship between the current density, J,and voltage, V has been derived by 
Soedergren et al. [26] based on electron diffusion and is given by: 


J =J, — 2 tanh (5) (exp 1 1) (15) 


Lp 


where the short-circuit current density, J,, is: 


Isc = STH (r + tanh (=) + sern] (16) 
Lp 


1-L2e2 cosh(+) 
Lp 


Here, ® is the light intensity; Lp is the diffusion length of the electrons; e the light absorption 
coefficient, and m is the ideality factor. The maximum voltage, V and current density, Jm, short circuit 
current density, Jsc, and open circuit voltage, Vc are extracted from the J-V curve and used to calculate the 
fill factor, FF, and the efficiency, 7 of the PSC by (17) and (18): 


— VmJm 
Bie IscVoc gy 
n(%) = Podst y 100 (18) 
in 
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2.3. Numerical simulation: method of line technique 

The MOL technique converts partial differential equations (PDEs) into a set of ordinary differential 
equations (ODEs) [16], [27]. The ODEs are obtained by discretizing the derivatives of spatial variables in 
(11) to (16) using finite difference approximation with applied boundary conditions. The parameters of the 
perovskites solar cell are listed in Table 1 [24], [26], [28]. The spatial variables in (9) to (14) are discretized 
using the centered difference approximation in (19) to (21) which offers a small tolerance for error and high 
accuracy. 


Ov; 1 ( +p; )( z = i)-¢ +p; dC =Q; ) 

5 = a [Piss = 2p; + pi- + PitPi+1)($i+1-Êi : PitPi-1) (Pi-Gi-1 | + y(nipi — 1) (19) 

i = = [Piss — 2n; +n; — TiTa Ai E a e e] +y(n;p; — 1) (20) 

ap; 1 1 

oe T p Piss — 26: + Qiz] — z u — p:i) (21) 
These equations can be expressed in general matrix form: 

M u = Au + f(x, t) (22) 


In (19) and (21) is transformed into matrix form shown in (22) for i = 1,2,3 ...m. The solution 
vector denotes as u, the time-dependent mass matrix, M, (where 1 denotes the time derivative and 0 denotes 
no time derivative), A is the discretization of the differentiation matrix in (9) to (11) with the boundary 
equations from (12) and (13), while f represents the generation and recombination of the charge density of 
holes and electrons in (10) and (11) and the net of charge densities in (9). The ODEs (19) to (21) in the 
matrix form (22) is then solved by using the implicit ODE solver ‘ode/5s’ in time forward via MATLAB 
with the resolution of the system, m = 101. Table 1 shows the PSC parameters for the simulation. 


Table 1. The parameters of the PSC [24], [27], [28] 


Parameter Symbol Value Unit 
Permittivity E 10008, A?stkg “m~ 
Diffusivity of holes Dp 2.5 x 107° m?s! 
Diffusivity of electrons Dn 2.5 x 107°? ms-* 
Boltzmann constant k 1.38 x10°23)— m@kgs~?K7+ 
Elementary charge q 1.603 x 10719 C 
Thickness L 0.5 x 1076 m 
Radiative recombination coefficient a 1.1 x 1078 m-3s71 
Light intensity © 1.5 x 107+ mst 
Diffusion length Lp 1.4 x 10% m 
Light absorption coefficient € 1.2 x 107 m+ 
Ideality factor m 1.3 


The dimensionless initial conditions are fitted from [24] resulting in (23) for holes and (24) for electrons: 
B(%,0) = 315.9x7 + 142.2x° — 71.48x> — 35.19x4 + 6.079x3 + 2.251x? — 0.06481x + 1.05 (23) 


ACZ, 0) = 175x® — 62.77x7 — 70.01x° + 3.829x° + 7.253x* — 1.403x3 — 0.2368x? — 0.006944x + 1.064(24) 


3. RESULTS AND DISCUSSION 

The result and discussion for the charge densities against thickness and the effect of the thickness 
variation on the efficiency of the PSC are discussed in this section. Figure 2 discovers the comparison of the 
simulation results (solid lines) and experimental data (circles) [24] for dimensionless electrons (blue color) 
and holes (red color) densities for L = 0.5 um. 

Figure 2 shows the comparison of holes (p) and electrons (n) densities between the experimental 
data [24] and simulation results obtained for 0.5 um cell thickness in dimensionless form. The simulation 
compares reasonably well to the experimental data. At X = —0.5 (the interface boundary between the hole 
transport layer (HTL) and perovskite layer), the electron density is increased, and the hole's density is 
decreased. At the interface of perovskite and HTL ¥ = —0.5, the charge separation occurs which allows only 
holes to get into the HTL layer while the electron remains in the perovskite layer. Meanwhile, it is vice versa 
for the density profiles % = 0.5 (the interface boundary between the electrons transport layer (ETL) and 
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perovskite layer). The charge accumulation model can be shown at the interface of the perovskite layer with 
the electrons transport layer (ETL) at X = 0.5 and the holes transport layer (HTL) at X = —0.5. 


3r 3 
—-—-—- p simulation at L=0.5um 
Biss n simulation atL=0.5um 
25+ O  nreference at L=0.5um [25]| Ọ2 5 
O preference at L=0.5um [25] | į 
i 
2 ? 419 
= ọ p a 
tt t+ 
af 1.5 R P Jisan 
i X g = 
i oe i 0” 
4 se IG OO ETAL, 1 
ao O~., 
“6 x, 
> 
a, 
0.5 -F 9 0.5 
0 + 0 


Figure 2. The comparison of the simulation results (solid lines) and experimental data (circles) [24] for 
dimensionless electrons (blue color) and holes (red color) densities for L = 0.5 um 


3.1. The effect of thickness variation of perovskite layer upon cell performance 
The perovskite layer thickness is varied from 0.15 wm to 0.45 um as referred to [29]. The value of 
the dimensionless parameters based on thickness is listed in Table 2. 


Table 2. The dimensionless parameters of the PSC 


Thickness, L (um) y À 
0.15 1.1 x 1073 1.004 x 1077 
0.25 3.0 x 107? 3.62 x 107? 
0.35 5.8 x 107 1.84 x 107? 
0.45 9.6 x 107? 1.12 x 107? 
0.5 1.18 x 107? 9 x 1073 


where A and y are the ratios of the Debye length to the width of the perovskite layer and the ratio of the 
typical charge diffusivity in bulk density to the typical diffusivity of the charges in the perovskite. Based on 
Table 2, the thickness of 0.15 um has the lowest value of y , that is 1.1 x 10-3 and increases to 9.6 x 107? 
as the thickness increases to 0.45 um. It shows that the recombination factor across the diffusion length 
slightly increases as the thickness increases in (10) and (11). At 0.15 um, the value of A is 1.0004 x 1071 
and reduces to 1.12 x 107? when thickness increases to 0.45um demonstrates that the ratio of Debye length 
(charge separation length) to the diffusion length decreases as thickness increases in (9). It allows more 
charges to diffuse toward the interface boundary (HTL/perovskite and perovskite/ETL). 

Figure 3 shows the variation of thicknesses are 0.25 um, 0.35 um, 0.45 um and 0.5 um which were 
scaled to dimensionless thickness, X = —0.5 to X = 0.5 by the scaling of (7). The results indicate that the 
higher the thickness, the electron and hole densities are slightly lower. For instance, a cell thickness of 0.15 
um has the highest electron and hole density value compared to a cell thickness of 0.45um which has the 
lowest electron and hole density value because of the increasing recombination process that can be related to 
y in scaling from (14). A higher thickness of perovskites affects the ratio of recombination rates, thus 
reducing the charge densities of PSCs [30]. The reduction of the charge densities occurs when the 
perovskite’s thickness is higher (the relation can be seen in (2) and (3)) and produces lower current density. 

Figure 4 discovers the effect of varied thickness on the concentration of the charge densities as time 
increases. The charge densities of the holes and electrons are more concentrated at L = 0.25 um refer Figure 4(a) 
compared to the charge densities at L = 0.5 um as shown in Figure 4(d). At the interface of the perovskite and 
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holes transport layer, the density of the holes becomes steeper as the thickness from L = 0.25 um increases to 
L = 0.5 um (refer Figure 4(a)-(d)) due to the effect of the Debye layer [30]. The situation at the interface of the 
perovskites and electrons transport layer also can be seen that the electrons densities become steeper from 
thicknesses L = 0.25 um to L = 0.5 um due to decreasing ratio of Debye layer. The permittivity gives an effect 
to the potential difference, @ in (2) and (3) for holes and electrons and results in a different level of current density. 
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Figure 3. The dimensionless electron and hole densities for various thicknesses 
(0.15 um, 0.25 um, 0.35 um, and 0.45 um) at time t = 5 x 1076s 
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Figure 4. The dimensional holes and electron densities for: (a) L = 0.25 um, (b) L = 0.35 um, 
(c) L = 0.45 um, and (d) L = 0.5 um 
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Figure 5 illustrates the current density, J, and voltage, V curve for 0.15 um, 0.25 um, 0.35 um, 0.45 um 
and 0.5 um. The maximum voltage, V„ and current density, Jm, short circuit current density, J/,. and open 
circuit voltage, Vç are extracted and tabulated in Table 3 which demonstrates the efficiency of different 
thicknesses. The value of the efficiency of the PSC is the lowest at thickness 0.15 um, which is 19.19% 
because the electron and holes move toward respective layers with too short distances with high charge 
densities. The efficiency started to increase from 19.19% to 19.77% as the perovskite thickness increased 
from 0.15 um to 0.25 um, due to high charge densities of holes and electrons and reduced the recombination 
and increasing ratio of Debye layer. However, the efficiency of the perovskites started to reduce from 
19.33% to 19.08% at a thickness of 0.35 um to 0.5 um, because the holes and electrons started to recombine 
while traveling toward selective layers respectively and reducing the ratio of Debye layer. As a result, the 
charge densities of holes and electrons started to decrease shown in Figure 4. Based on the efficiency result 
on different thicknesses in Table 3, showing that the efficient thickness of perovskite ranges from 0.25 um. 
Hence, thickness plays an important role in affecting the performance of the solar cell. 


9 


Current density Am “ 


50 F l| 4 
| 
il 
4\|| 


0 0.2 0.4 0.6 0.8 1 1.2 
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Figure 5. The J-V Characteristic curve for different thicknesses 


Table 3. The calculation of efficiency of the perovskites solar cell based on different thicknesses 
Thickness, L (um) Js (AM?) VeV) Jm (Am™) WU) FF 1% 


0.15 223.92 1.01 213.175 0.9 0.8483 19.19 
0.25 237.59 0.99 219.663 0.9 0.8405 19.77 
0.35 239.87 0.98 214.772 0.9 0.8223 19.33 
0.45 240.251 0.97 238.607 0.8 0.8191 19.09 
0.50 240.296 0.97 238.469 0.8 0.8185 19.08 


4. CONCLUSION 

A drift-diffusion model of a one-dimensional solar cell in the perovskite layer is presented and is 
solved numerically. The charge densities of holes and electrons' behavior across the thickness are shown. The 
smaller the thickness of the cell the lower the cell efficiency. The best efficiency obtained from the result is 
19.77% at a thickness of 0.25 um which is proven that the optimum thickness to achieve maximum efficiency, 
but the thickness of perovskite is still acceptable until 0.5 um. The model would be useful as a guideline for the 
thickness of cells that predicts the perovskite cell performance. However, the simulation can be improved by 
investigating the behavior of charge density across the three layers of the solar cells. Therefore, the folding 
technique is suggested to improve the numerical scheme of holes and electrons for all phases for future work. 
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